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Abstract 

In a semiconductor superlattice with long scattering times, damping of Bloch oscillations due 
to scattering is so small that convective nonlinearities may compensate it and Bloch oscillations 
persist even in the hydrodynamic regime. In this case, numerical solutions show that there are 
stable Bloch oscillations confined to a region near the collector with inhomogeneous field, charge, 
current density and energy density profiles. These Bloch oscillations disappear when damping due 
to inelastic collisions becomes sufficiently strong. 
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I. INTRODUCTION 



Bloch oscillations (BOs) are coherent oscillations of the position of electrons inside en- 
ergy bands of a crystal under an applied constant electric field —F. Their frequency is 
ub = eFl/h (/ lattice constant), and therefore it can be tuned by an applied voltage. BOs 
were predicted by Zener in 1934 as an immediate consequence of the Bloch theorem 
but they were not observed for almost sixty years because scattering of the electrons with 
phonons, impurities, etc. damp them very rapidly into oblivion. To observe BOs, their pe- 
riod has to be shorter than the scattering time r, and therefore the applied field has to 
surpass the value /?-/(e/r), which is too large for most natural materials, in which / is of 
Angstrom size. In 1970, Esaki and Tsu suggested to create an artificial crystal, which they 
called superlattice (SL), by growing many identical periods comprising a number of layers of 
two different semiconductors with similar lattice constants The period of the resulting 
one-dimensional crystal may be much larger, say about 10 nm, and this gives reasonable 
electric fields of about 10 kV/cm, which are within the range of experimental observation. 
Damped Bloch oscillations were first observed in 1992 in semiconductor SLs whose initial 
state was prepared optically Besides their interest for theoretical physics, BOs have 
attracted the attention of many physicists and engineers because of their potential for de- 
signing infrared detectors, emitters or lasers which can be tuned in the THz frequency range 
simply by varying the applied electric field j^]. However no electrically driven devices based 
on BOs have been realized. Another application is based on the fact that BOs give rise 
to a resonance peak in the absorption coefficient under dc+ac bias and a positive gain at 
THz frequencies jsj. The latter has been observed in quantum cascade laser structures jo]. 
These applications are severely limited by scattering which rapidly damps BOs and, for a 
dc voltage biased SL, favors the formation of electric field domains (EFDs) whose dynamics 

nn 

yields self-sustained oscillations of lower frequency (GHz) [7|, |8| (a phenomenon similar to 
the Gunn effect in bulk GaAs [9]). EFD formation may also preclude THz gain in simple 
dc+ac driven SL which is typically calculated assuming spatially uniform solutions of drift- 



diffusion or Boltzmann type equations 
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14| . This assumption has not been tested by 



solving space-dependent equations with appropriate boundary conditions or by experiments 
in semiconductor superlattices. An interesting idea for efficient terahertz harmonics gen- 
eration is to excite relaxation oscillations in the superlattice by incident radiation from a 
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waveguide |15 |. 

To understand the role of EFD formation in the observation of BOs or THz Bloch gain, our 
starting point should be a model in which BOs and EFDs are both possible solutions of the 
governing equations. One simple possibility is to use a self-consistent version of the Ktitorov, 
Simin and Sindalovskii (KSS) Boltzmann equation 151 with Bhatnagar-Gross-Krook (BGK) 
collision terms |16|, 117[. The characteristic equations of this kinetic equation exhibit BOs 
as solutions, whereas there exist a hydrodynamic regime for large applied electric fields 
;hat yields Gunn-type oscillations of the current as solutions of a drift- diffusion equation 
ItI - However the KSS model does not have a hydrodynamic regime of the kinetic equation 
in which both BOs and EFDs are possible solutions. Why is the coexistence of BOs and 
EFDs not possible in the KSS model? Firstly, the BOs and the hydrodynamic regime 
correspond to widely separated time scales. In the slower hydrodynamic regime, BOs have 
already disappeared due to scattering and we cannot study simultaneously BOs and the 
EFDs appearing in the hydrodynamic regime. This problem could be reduced if we consider 
materials with long-lived BOs corresponding to almost elastic collisions, as will be discussed 
in more detail later. The second problem is that, unlike the original BGK local equilibrium 



181], the local equilibria used in 16| and in |tl7|] depend only on the 2D electron density 
n, which remains approximately constant during a BO. Since hydrodynamic regimes are 



perturbations of local equihbria, the KSS model as modified in [17| can only provide slowly 
varying drift-diffusion equations 19| for the electron density and the electric field —F that 
cannot contain BOs among their solutions. This second problem can be solved if the local 
equilibrium in the Boltzmann-BGK kinetic theory depends on electron density, electron 
current density and mean energy and the collision term preserves charge but dissipates 
momentum and energy [2^. The most important property of the proposed model is that 
it allows the local equilibrium distribution to oscillate at the Bloch frequency, which is the 
crucial feature (missing in the KSS kinetic equation) if we want to derive a hydrodynamic 
regime that allows BOs. 
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II. MODEL 



The model equations are 



dtf + v{k) dj + eFh-'duf = Q[f] = -y{f - /^), (1) 
ed^F = er\n-ND), (2) 



^ ^ukl+P cos kl 



f''(k;n,Jn,E)=n——- , (3) 



n 



— f(x,k,t)dk = — fdk. (4) 



f{x,k,t)dk = — 

t/I J-n/l 

Here n, N^, e, — e < 0, m*, u, and —F are the 2D electron density, the 2D doping density, 
the permittivity, the electron charge, the effective mass of the electron, the constant collision 
frequency and the electric field, respectively. v{k) = Al sm{kl) / {2h) is the group velocity 
corresponding to the miniband tight binding dispersion relation £{k) = A(l — cos kl)/2. For 
the sake of simplicity, we have assumed a Boltzmann local equilibrium (j3]), but it is easy 
to replace it by the Fermi-Dirac local distribution in the degenerate case. The distribution 
functions / and have the same units as n and are 27r//-periodic in k (the function ukl 
in (E]) is extended periodically outside — vr < kl < n). A quantum version of (II]) can be 
obtained as indicated in Ref. |8|. 

The dimensionless multipliers /3(x,t) and u{x,t) depend on J„=e ^^^^^viji) f dk/ {27^) 
(electron current density) and on E=l J^^|,JA/2 — £{k)]fdk/{27m) (mean energy). They 
are found by solving 
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■k/I 




v{k) f 


dk 


-TX/I 












l-n/l V 2" " 





(5) 



Similarly to BGK collision models in rapid granular flows 2l|], the restitution coeffi- 
cients aj and ae take values on the interval [0, 1] and measure the dissipation due to 
collisions in current density and energy, respectively. In fact, the collision operator 
satisfies f^^J^^Qlf] dk = (charge continuity), e f^^J^^v{k) Q[f]dk/{2TT) = —uajJn, and 
/ jZi^ii^f^ - S{k)]Q[f]dk/{2nn) = -ua^iE - Eq). Obviously for a^j = the collisions 
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conserve energy and momentum (elastic limit). To simplify matters, we shall assume that 
aj and are constant. Eq is the mean energy at the lattice temperature of the global equi- 
librium which will be reached in the absence of bias and contact with external reservoirs. 
Equivalent results are obtained if we define the mean energy as the average of £{k), which is 
equal to A/2 — E', but Eq. ^ leads to a simpler relation between energy density and lattice 
temperature. At the lattice temperature, Tq = A/(2A;b/3o), m = 0, = i?o, and ([5]) yields 
2Eo/ A = Ji(/3o)//o(/3o)) where Is{x), s = 0,1, are modified Bessel functions. 



III. HYDRODYNAMIC EQUATIONS 



In the hyperbolic limit in which the collision and Bloch frequencies are comparable and 
dominate all other terms in ([1]), it is possible to derive closed equations for nondimensional 
n, F and A (the complex envelope of the BO solution) 22|, 



fi = A{x,t)e-'' + fi,s{x,t) 



F{x, s) ds 



(6) 



{6 is the rapidly varying phase of the BO, fi^s = 0{5) is written below), provided the 
collisions are almost elastic. The small dimensionless parameter 5 = e^Nj:,l/\/ {2efi^v'^) is 
the ratio between the scattering time and the dielectric relaxation time and the restitution 
coefficients are assumed to scale with it, = 57ej- A similar double limit of vanishing 
Knudsen number (equivalent to 5 — )■ 0) and almost elastic collisions has been used to derive 



23II24. In ([6]), fj{x,e,t;6) are the Fourier 



hydrodynamic equations for granular gases 
coefficients of f{x, k, t; 6) = J^'jL^oo fj^^^^ A — nE — iJn. The nondimensional equations 
are 



dF 
'dt 



+ 



TP TP , ^ 1 •'2,0 

'ypEnnr H T^lm ■ 

' " 2 dx l + 2iF 



rB{0) 
J2,0 

l + 2iF 



OF 
dx 
dA 
'dt 




- FRe hs 



+ 57elm/i5' 



(7) 
(8) 

(9) 
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h. = 



fi^su dlmfi 



5n 



n 



f 



l,Su 



dx 

5'jenEo{5'yj - iF) 



(10) 



1,5 



nEs - iJn,s 
6 



[7enEo(57j - iF) 



+ 



F^ + .7, 

- iF)Rehs - {F + i6-fe)'^T^hs 

F + d ( ^ f,T 

— n — Re ■ — 

2 dx \ l + i2F 



+ — — Im- ' 



dx \1 + i2F 
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We have defined the nondimensional variables / = f /No, n = u/Nd, E = 2E/A, Jn = 
J/[Jn], X = x/[x], ... (where [y] are the units in Table H]) and omitted tildes over variables. 
The dimensionless multipliers /3 and u in are functions of the rapidly varying BO phase 
9 due to and therefore, we can expand in ([3]) in powers of 6, ~ -\- Sf^^^\ 
The /^("") (m = 1,2) are now 27r-periodic functions of and k. Then we have the Fourier 
coefficients 



f 

J 1 



B{0) 



-ijk—imO 



— TT — TT 



dkde 



(12) 



in which we set /i = Ae""^^ ignoring 0{5) terms in ([6]). 

To derive ©-Q, we start from the equations for the moments fj which can be obtained 
from ([1]) by integration over k [2(|: 



— Im— — 

ot ox 



0, 



iF]h = 5 



' f p 7e + 7i f le - 7j p 1 ^ r f f 



(13) 

(14) 



where /o = n, fi = nE — iJn, and there are similar equations for higher moments. From f[T3|) 
and the Poisson equation dF/dx = n — 1, we find Ampere's law for F: dF/dt = J{t) — J„, 
where J{t) is the total current density. We shall assume that the second moment /2 is a 
known function of /o and /i, /2 = 5'(/o, fi)- Then we find equations for n, F and A in (E]) by 
a method of nonlinear multiple scales [25| with time scales 9 and t. To obtain the function 



(7, we carry out a Chapman- Enskog expansion 



13, 



22| for ([T]) (in dimensionless units) with 
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a time derivative given hj F df/d9 + Sdf /dt, according to (Q. The distribution function 
is supposed to be periodic in k and in 6. This procedure gives approximate formulas for 



g = f2 from which ([7])-© are obtained 22|. 



The hydrodynamic equations ([7])- OH]) have the spatially uniform solutions, n = 1, 
J = 6'jeEQnF/ {6'^'~fe'~fj + F"^ ) (i n dimensional units this gives the well-known temperature 



dependent drift velocity 16| 26|), and A = A^e'^'^^^"'^^^^'^. Inserting the latter formula in (|6]), 
we see that this corresponds to a damped BO whose amplitude relaxes to 0. Even when we 
manage to prepare the initial state with a coherent BO of complex amplitude Aq, ignoring 
space dependence will lead to disappearance of the BOs after a relaxation time 2/(7e + 7j). 
Stabilization of the BOs may be caused only by the spatially dependent second term on the 
right hand side of OH])- 

TABLE I: Hyperbolic scaling and nondimensionalization with = 10^^ Hz. 
/, n F £, E v{k) Jn x k t 6 

el 2 2h 2h eWB I e^NolA 2eK^u^ 

IQio kV 10''™ lO^A „^ ^ „„ 

cnP' cm s cnP" nm 

4.048 130 8 6.15 7.88 116 0.2 1.88 0.0053 



IV. RESULTS 



We now solve numerically the hydrodynamic equations with the boundary conditions [8| 



a 



dF 
'dt 



+ (roF 



J, 



x=0 



OF 
'dt 



+ aiuF 



J, 



x=L 



1 

— I F{x, t) dx = 0, 
L Jo 



ox 



at X = 0. 



(15) 

(16) 
(17) 



We obtain similar numerical results with A = a.t x = 0. Here L = Nl/[x] and (p = 
eV/ {huN) are the dimensionless SL length and average field (proportional to the applied 
voltage V), respectively. We have used contact conductivities cro,i = 12.1 (fim)~^ which 
yield dimensionless conductivities ctq,! = 0.2 (conductivity units are [a] = e^Ni:,M/{2h^i')). 
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Initially, F{x,0) = and A{x,0) = Aq (constant). The latter condition means that we 
have prepared the SL in an initial state having a coherent BO with complex amplitude Aq. 
Whether this can be achieved by optical means as in js] remains to be seen. 

We solve ([T])-® with the parameter values indicated in Table |I] (which are similar to those 



in Ref. 



27j ) and different values of the restitution coefficients. We start with ctg = aj = 0.01 



so that z/cte = uaj = 10^^ Hz. 28(] The 3D doping density N^d = 8 x 10^^ cm ^ gives 
Nd = NsdI = 4.048 X 10^° cm'^ as in Table [D and e = 12.85 eo- We find 6 ^ 0.0053 
and 7ej = dej/S = 1.8781. We consider a 50-period (A^ = 50) dc voltage biased SL 
with lattice temperature 300 K. For V = 0.2 V (therefore = 0.06 I29I), we observe that 
|y4(x, t)| first diminishes uniformly from Aq = 0.153 to almost zero after a relaxation time 
2/(7e + 1j) ~ 0.53 (about 1 ps). Later a small pulse is formed at about x = L/4 which 
subsequently extends to the remaining part of the sample and it grows nearer to its end. 
The BOs are confined to the second half of the sample that is closer to x = L and are zero 
in the first half of the sample closer to x = 0. Thus the profile of |A| has a compact support 
with a maximum near x = L. \A{x, t)\ is close to a periodic oscillation in time: small pulses 
are formed at the left of its support, climb up towards the maximum of the pulse which 
then diminishes and the same behavior repeats itself. Figure [T]^a) shows four snapshots of 
|y4(x, t)| illustrating this behavior which can also be observed in two movies 30|. The field 
profile depicted in Fig. [IJb) is almost stationary. The mean energy and electron current 
densities during one BO can be reconstructed by means of ([6]). We show them for = 
in Fig. [U^c). At the two different SL locations marked in Fig. [H^b), the graphs of Jn versus 
time are shown in Fig. Wl^d). This figure and two additional movies showing the evolution of 
the Jn and 1 — E profiles in the supplementary material illustrate that the Bloch frequencies 
depend strongly on space and are higher near the collector where the field is larger. 

For the scattering times reported in Ref. [27|, the restitution coefficients are = 0.09 
and aj = 0.29, but the BO amplitude becomes zero everywhere after a short relaxation time. 
BOs also disappear for = 0.01 and aj = 0.29/9 ~ 0.032, and there is a smaller critical 
value of aj (for fixed ae) below which BOs can be sustained. They also persist for aj = 0.01 
and ae = 0.09/29 which keep the same ratio ajja^ = 29/9 ~ 3.22 as in Ref. j27|. There 
is a critical curve in the plane of restitution coefficients such that, for {■y^ + 7j)/2 > 7crit 
(7crit ~ 2.5 for 6 = 0.0053), BOs disappear after a relaxation time but they persist for 
smaller values of (7e + 7j)- 
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FIG. 1: (a) Modulus of the BO complex amplitude vs space at times ti = 7, t2 = 9, ^3 = 11, 
t4 = 13. (b) Stationary field profile, (c) Profiles of the nondimensional mean energy density E 
and nondimensional electron current density J„ for 9 = 0. (d) Current density at the two different 
points marked by (1) and (2) in (b) during BOs. Clearly, the frequency at point (2) is larger than 
at (1). To transform the magnitudes in this figu9e to dimensional units, use Table HI [A] = ANd/2. 



In summary, we have analyzed the Boltzmann-BGK-Poisson equations with local equilib- 
rium depending on the electron density, current density and energy density in the hyperbolic 
limit in which the BO period is much shorter than the dielectric relaxation time and colli- 
sions are almost elastic. In the long-time scale, there is a hydrodynamic regime described 
by coupled equations for the electric field, the electron density and the BO complex ampli- 
tude. When the restitution coefficients (equivalently the inverse of the scattering times) are 
sufficiently small and the initial state has been prepared so that there is a nonzero Bloch 
oscillation, there are stable spatially inhomogeneous profiles of current and energy densities 
displaying BOs confined to a fraction of the SL extent. It remains to investigate whether 
confined Bloch oscillations may also be important for terahertz harmonic generation with 
underlying inhomogeneous charge and electric field profiles. 
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